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Abstract 



We explore theoretically the principles that govern photon emission from 
single-molecule conductors carrying electric currents between metallic con- 
tacts. The molecule and contacts are represented by a generic tight-binding 
model. The electric current is calculated using Landauer theory and the pho- 
ton emission rate is obtained using Fermi's golden rule. The bias-dependence 
of the electronic structure of the molecular wire is included in the theory in 
a simple way. Conditions under which significant photon emission should oc- 
cur are identified and photon spectra are calculated. We predict the photon 
emission rate to be more sensitive than the electric current to coupling asym- 
metries between the molecule and contacts. This has important implications 
for the design and interpretation of STM experiments searching for electro- 
luminescence from individual molecules. We discuss how electroluminescence 
may be used to measure important characteristics of the electronic structure 
of molecular wires such as the HOMO-LUMO gap and location of the Fermi 
level of the contacts relative to the HOMO and LUMO. The feasibility of ob- 
serving photon emission from Au/benzene-dithiolate molecular wires is also 
discussed. 

PACS: 78.67.Lt, 78.60.Fi, 73.63.Rt 
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I. INTRODUCTION 



A molecular wire is a single molecule (or a few molecules) that forms an electrically 
conducting bridge between a pair of metallic nano- contacts. In recent years molecular wires 
have been realized in the laboratory and their electronic transport properties have been 
measured 1-7 . There have also been many theoretical advances in modeling such systems 2 ' 8-31 
and the study of transport in molecular wires continues to be an active area of experimental 
and theoretical research 32 . It has been found that the current flowing through a molecular 
wire depends strongly on the electronic structure of the molecule as well as the geometry 
of the molecule and contacts and the chemical bonding between them 8-15 ' 19 ' 22-25 ' 27-31 . In 
a simplified picture, when a potential bias is applied to the contacts, the electrochemical 
potentials of the contacts separate and molecular orbitals located in the window of energy 
between the two electrochemical potentials mediate electron flow from one contact to the 
other. The current rises sharply each time this window of energy expands to include an 
additional molecular orbital. 

Much work has been done modeling and attempting to understand I-V characteristics 
that have been obtained experimentally 2 ' 14 ' 17 ' 19 ' 21 ' 24 ' 26 ' 30 . No theoretical studies have been 
reported, however, of another potentially important property of self-assembled molecular 
wires, namely, photon emission (electroluminescence) from a molecular wire carrying an 
electric current due to electronic transitions between molecular orbitals. When a bias volt- 
age is applied to the contacts, the molecular wire moves out of equilibrium, with a flux of 
electrons passing through it. Electrons enter the molecule from the contact with the higher 
electrochemical potential, and drain into the other contact. Since the electrochemical poten- 
tials of the contacts are no longer the same, one can no longer assume that the orbitals of the 
molecule will be filled up to an energy equal to a common Fermi energy of the contacts. The 
passage of electrons through the molecule may lead to partial occupation of various different 
molecular orbitals that lie within the electrochemical potential window of the contacts. It 
is possible that transitions from one partially occupied molecular orbital to another of lower 
energy occur, resulting in photon emission. 

Although definitive experimental evidence of this effect has not yet been reported, re- 
lated phenomena have been observed in scanning tunneling microscopy (STM) experiments: 
Systems with an STM tip over a clean metallic surface are known to emit photons due to 
the decay of plasmons 33 . Recent STM experiments on molecular monolayers adsorbed on 
metal substrates have suggested that a different, molecule-dependent photon emission pro- 
cess may also occur: Poirier has reported molecule-dependent photon emission in an STM 
experiment involving reduced and oxidized alkanethiol monolayers adsorbed on Au(lll) 34 . 
In an STM experiment involving monolayer films of Cqq fullerenes on Au(llO) surfaces, 
Berndt et al. observed enhanced photon emission when the STM tip was placed above an 
individual molecule 35 . Smolyaninov has reported photon emission at energies similar to cop- 
per phthalocyanine (CuPc) transition energies in STM experiments involving CuPc adsorbed 
on Au 36 . However, in very recent STM experiments by Hoffmann et al. on hexa-tert-butyl- 
decacyclene on noble-metal surfaces 37 , it was found that while the presence of the molecules 
modulated the photon emission, the observed emission appeared to be plasmon-mediated 
rather than molecular in character. Thus, the nature of photon emission mechanisms for 
STM-molecular monolayer systems is yet to be adequately understood. 
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In all of these STM experiments, a bias voltage was applied between the STM tip and 
metallic surface, with a molecular monolayer interposed between these two contacts as in 
molecular wire systems. However, there are important differences between these STM molec- 
ular monolayer systems and self-assembled molecular wires: The coupling between the STM 
tip and the molecule is usually very small compared to the coupling between the molecule 
and the substrate, whereas for self-assembled molecular wires that are bonded chemically 
to both metal contacts, the couplings of the molecule to the source and drain may be com- 
parable in strength. Also, the STM experiments involve planar metallic substrates, which 
are quite different from the nanoscopic metallic contacts of some self-assembled molecular 
wires. It is unclear how these differences would affect photon emission, so it is of interest to 
consider photon emission not only from STM-molecular monolayer systems but also from 
more symmetrically coupled self-assembled molecular wires. 

Electroluminescence from molecular wires is a potentially important effect, not only be- 
cause of its intrinsic fundamental interest, but also as a novel experimental probe of molecular 
nano-electronic devices. For example, as we will show below, by observing the photon emis- 
sion spectrum of a molecular wire, important information about its the electronic structure, 
including the energies of the molecular orbitals, their locations relative to the Fermi level, 
and their occupations as a function of bias voltage, may in principle be obtained. There is 
currently disagreement between different theoretical models of 1,4 benzene-dithiolate (BDT) 
attached to gold contacts 14 ' 19 ' 21 ' 25 ' 28 ' 30 regarding the energies of the highest occupied molec- 
ular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), relative to the 
Fermi level of the contacts. If the photon emission spectra of such systems were modelled, 
and compared with experiment, important new insights into the most appropriate models 
for molecular wires could be obtained. 

In this article, we present and solve a generic model for photon emission from a molecular 
wire. As this is the first theoretical study of this effect, our purpose will be to develop a 
basic qualitative picture of the underlying mechanism, and to examine how adjustment of 
the various model parameters, to reflect different experimental situations, should affect pho- 
ton emission. We use a one-dimensional tight-binding model that attempts to capture the 
important physics involved in a way that is simple and intuitively reasonable, and is qual- 
itatively consistent with what is already known about molecular wire electronic structure 
and transport. Another objective of this work is to provide some guidance to experimen- 
talists as to which types of systems to study in order to observe photon emission. We find 
that, under certain circumstances, the model predicts that photon emission should occur, 
due to molecular orbital transitions. We show how molecular orbital occupations behave 
as a function of applied bias voltage, and how varying the model parameters changes this 
behaviour. Fermi's Golden Rule is used to calculate photon emission spectra. We show how 
the emission depends on molecular orbital occupations, and the locations of the molecular 
orbitals relative to the Fermi levels of the contacts. 

Our model and theoretical approach are described in Section II. The results of our cal- 
culations and their interpretation are presented in Section III. We conclude in Section IV by 
summarizing our results, and discussing their implications for possible experiments directed 
at observing photon emission from BDT attached to gold contacts. 
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II. THE MODEL 



Each metal contact will be modelled as a one-dimensional tight-binding chain. The 
molecular wire is modelled as a pair of atoms placed next to the origin, forming a bridge 
between the two contacts (see Fig. 1). The model Hamiltonian of this system is 

H = ^e L \n){n\ + (3{\n - l)(n\ + \n)(n - 1|) 

n<0 

+ ^2e R \n){n\+ P(\n + l)(n\ + \n)(n + l\) 

n>0 

+e a \a){a\ + e b \b)(b\ +/3_i, (| - l)(a| + |«><— lj) 

+P 1 , b (\l)(b\ + \b){l\)+Pa,M(b\ + \b){a\) (1) 

where ep and ep are the site energies of the left (source) and right (drain) contacts, and 
e a and e b are the site energies of atoms a and b of the molecule, which depend on the bias 
voltage applied to the wire and the nature of the metal contacts, as well as the identities 
of atoms a and b. j3, /3-i, a , Pi,b and f3 a>b are the hopping amplitudes between atoms in the 
contacts, between atom a of the molecule and the left contact, between atom b and the right 
contact, and between atom a and atom b, respectively. \n) represents the orbital of an atom 
in one of the contacts, and \a) and \b) represent the atomic orbitals of the molecule. We 
consider the orbitals of different atoms to be orthogonal. However, this model may easily be 
extended to systems where this is not the case 16 . We take the the electrochemical potentials 
of the source and drain contacts to be us = Ep + eVbi as /2 and fiD = Ep — eVuas/2 where 
Vbias is the bias voltage applied between them and Ep is their common Fermi level at zero 
applied bias. The applied bias also affects the site energies e L and e R of the contacts so 
that e L = e contacts + eV bias /2 and 

^contacts - eV bias /2, where 

^contacts is the zero bias 

site energy. Our treatment of the effect of the bias voltage on site energies e a and e b of the 
molecule itself is described at the end of Section II. 

Electrons exist in the form of Bloch waves in the contacts, and undergo reflection or 
transmission when they encounter the two-atom molecule. Their wavefunctions are of the 
form 

|^) =^(e ifc " d + re- lfend )|n) + ^te lfc '" d |n) + Ca|a) + c 6 |6) (2) 

where d is the lattice spacing, and t and r are the transmission and reflection coefficients. 
The electrons have eigenenergies of the form E = q + 2(3cos(kd). This equation holds 
for both €i = €l and = €r, so when an electron with initial wavevector k undergoes 
transmission, its wavevector changes (to k') due to the difference between ep and ep. For 
certain k, there are no real solutions for k! . In these cases, k' becomes complex, and the 
transmitted Bloch wave is evanescent. By applying ( — 1|, (a\, (b\, and (1| to H\ip), analytic 
expressions for the transmission and reflection coefficients are obtained: 

2ip_ l!a p ljb p a>b (3sin(kd) 



r = 



\P(E - e a ) - (3\ a e-^][(3{E - e b ) - - (3^ 

[(3{E - e a ) - 0l^d**\\P(E ~ ^) ~ ^ ik ' d \ ~ PIP 2 



(4) 
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The coefficients c a and q> for the atomic orbitals in the molecule are similarly obtained, in 
terms of the transmission and reflection coefficients: 



1(1 + (3t 

C a = —5 , Cb = "5— (5) 

P-l,a Pl,b 



The transmission probability is given by 



, 2 v(k') . , 2 sin(k'd) 



T = 1*1 = 1*1 — TwV ( 6 ) 
v{K) sin{ka) 

where v(k) is the velocity of the incoming wave and v(k') is the velocity of the transmitted 
wave. In our model, T depends on the energy of the incoming electron, as well as on 
Vbias- The Vbias dependence arises from the effect of Vu as on the site energies of the contacts 
(described above) and of the molecule that will be discussed below. Using Landauer theory 38 , 
an expression for the current is obtained: 



2e r 



MS 

T(E,V bias )dE (7) 

fJ-D 



Here we assume the temperature to be K, so the Fermi functions are trivial, and the 
transmission probability is only integrated through the Fermi energy window. 

To calculate emission spectra for the molecular wire we use the expression for the spon- 
taneous emission rate of a system emitting photons into empty space, using Fermi's Golden 
Rule 39 . The emission rate is given by 

Ae 2 u 3 |2 



3ftc 3 



l(^/|x|^)| 2 (8) 



where ipi an d represent initial and final states, and hu is their difference in energy. We 
consider emission only from the molecular sites a and b. The emission rate is therefore 
approximated by 

4e 2 u 3 

R = l^l c «i/ Ca ^ a l x l a ^ + c h c bA b \ x \ b )\ ( 9 ) 

where i and / label initial and final states. The overlap terms (a|x|6) and (6|x|a) are 
neglected since they should be small compared to (a|x|a) and (6|x|6). We approximate 
(a|x|a) and (6|x|6) by the locations of their atomic centers, so that (a|x|a) = — |, (b\x\b) = | 
(b being the molecular bond length). Thus, we have 

c^uj^b 2 

R(ki,u) = -^-\ C b,f C b,i - C* aJ C ati \ 2 (10) 

To calculate the emission rate as a function of photon energy, we must consider all electron 
states of the system, incoming from both the source and drain contacts. Since we assume 
the temperature to be K, all states up to the electrochemical potential of the appropriate 
contact are occupied, ipf must be initially unoccupied, and it must be of lower energy than 
ipi. Therefore, we consider transitions from occupied initial states that are incoming from 
the source, to final states, within the electrochemical potential window, that are incoming 
from the drain. After normalizing the wavefunctions and converting the sum over k-states 



5 



(and spin) into an integral over energy, an expression for the photon emission spectrum (for 
a given bias voltage) is obtained: 



In order to interpret our results physically, it is useful to relate the emission spectra to 
total occupations of the molecular orbitals, as a function of bias voltage. For this purpose, 
we project eigenstates of the system given by equation (2) onto the bonding and antibonding 
orbitals of the isolated molecule, whose coefficients will be labelled Co and Co* respectively. 
Their total occupations are calculated by summing over all occupied electron states (in- 
cluding spin) incoming from both contacts. This sum is converted into an integral, and an 
expression for the total bonding and antibonding orbital occupations (O, O*) is obtained: 



Changes in the occupations of the bonding and antibonding orbitals in response to changes 
of the applied bias voltage may result in some charging of the molecule. If this occurs, the 
charging causes an electrostatic shift of the molecular energy levels that in turn severely 
limits the actual charging that takes place 23 . In the present work we approximate the 
shift of the molecular levels in response to the applied bias by adjusting e a and equally 
so as to maintain the net charge that the molecule has at zero bias. The simplicity of 
this approximation is in keeping with the generic nature of the molecular wire model that 
we consider here. It yields behavior of the molecular levels with bias that is physically 
reasonable, and as will be seen in Section III B, remarkably similar to that obtained from 
ab initio calculations for some molecular wire systems. 

In the remainder of this article we will consider the case where the Fermi level of the leads 
at zero bias falls within the HOMO-LUMO gap of the molecule, as is typical of molecular 
wires such as the Au/benzene-dithiolate/Au system. Thus we will identify the bonding and 
antibonding orbitals O and O* with the HOMO and LUMO, respectively, of the isolated 
molecule. (The terms O and HOMO will henceforth be used interchangeably as will O* and 
LUMO.) Our model parameters will be chosen so that at zero bias the Fermi level of the 
leads lies between O and O*. As will be seen below, if the Fermi level is well above O and 
well below O* (but O and O* are within the electronic energy band of the leads) it follows 
from Eq. 12 that at zero bias O is almost completely filled with electrons and O* is almost 
empty, although O is never totally full and O* is never completely empty because of the level 
broadening that occurs due to hybridization between the states of the molecule and leads. 
In this way the HOMO and LUMO character of the O and O* orbitals is expressed within 
our model when the molecule couples to the leads. By choosing different values of the site 
energies e a and e& at zero bias, the energies and occupations of the O and O* orbitals can 
be varied. Thus the effects of differing amounts of charge transfer between the molecule and 
leads can also be studied within our model. 

We also note in passing that although our results are described here in terms of electron 
transport and optical transitions between occupied and empty electron states, an equivalent 
description can be given in terms of electron and hole transport and optical transitions. 




(11) 




(12) 
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III. RESULTS 



We consider the case where the energy bands of the contacts are partially filled {k F = 
In numerical results presented below the hopping parameters are chosen to have values 
typical of molecular wire systems (f3 = —5 eV, f3 ajb = —1.5 eV), resulting in a HOMO 
LUMO gap of 3 eV for the isolated molecule. 

A. Symmetric molecule-contact couplings; HOMO LUMO gap centred at Ep. 

In this section, The HOMO-LUMO gap of the molecule is chosen to be centred at the 
zero bias Fermi level of the contacts, and the molecule-contact couplings have equal values. 
First, we consider the case where these couplings have values similar to those in chemically 
bonded molecular wire systems (/3-i, a = Pi,b — — 1-0 eV). Fig. 2a shows how the source 
and drain electrochemical potentials fis and /i£> and the energies of the bonding (O) and 
antibonding (O*) orbitals behave as a bias is applied to the contacts. At zero bias, O is 
located at -1.5 eV and O* at +1.5 eV. Fig. 2b shows the electron transmission probability 
through the molecule at zero bias obtained from eq. 6. The bonding and antibonding 
orbitals provide channels for electron transmission, so there are peaks at electron energies of 
-1.5 eV and +1.5 eV. The transmission peaks are somewhat broadened due to the coupling 
of the molecule to the contacts, which causes the discrete molecular orbitals (O and O*) 
to hybridize with the continuum of states in the contacts. Returning to Fig. 2a, as Vb ias 
increases, the electrochemical potentials of the source (/X5) and drain separate. At 

Vbias — 3 V, /is moves above O* and hd moves below O. In this symmetric case, this does 
not result in any change in the energies of O and O*. Fig. 2c shows the electron occupation 
of the molecular wire, as projected onto the bonding (O) and antibonding (O*) orbitals of 
the isolated molecule. At zero bias, O is almost fully occupied, and O* is almost empty. 
As Vbias increases, O partially empties. \id moves below O in energy, so electrons from 
the drain no longer contribute to the filling of the bonding orbital, while electrons from 
the source continue to contribute, so O becomes half-filled. Similarly, electrons from the 
source start to fill O* as us approaches O* in energy, and finally O* becomes half-filled as 
well. As O and O* empty and fill with V&j as , the total occupation remains constant (two 
electrons), due to the symmetric placement of E F , at the centre of the HOMO LUMO gap. 
Therefore, there is no tendency for the molecule to charge, and the orbital energies remain 
constant (Fig. 2a). Notice that the filling and emptying of O and O*, as a function of Hias> 
is gradual. This is due to the fact that the continuous density of states associated with 
the broadened resonances in Fig. 2b enters the electrochemical potential window between 
Us and /ip gradually. Fig. 2d shows the orbital occupation for the case of weak molecule- 
contact coupling (/3_i, a = Pi,b = —0.2 eV). The weak coupling results in a more sharply 
peaked molecular density of states, so O and O* empty and fill much more abruptly. Fig. 
2e shows the resulting total photon emission rate for (i) /3_i )0 = f3i^ = —0.2 eV and (ii) 
P-i,a — Pi,b = —1-0 eV. Emission is strong when the energies of O and O* are inside the 
electrochemical potential window. In case (i), emission increases more rapidly with bias 
near Vu as = 3.0 V than in case (ii), because O and O* empty and fill more abruptly. Fig. 
2f shows how, within our model, the emission spectrum for case (ii) changes with bias. As 
expected, emission is peaked around the transition energy of the molecule (3 eV) at higher 
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bias. (Weaker molecule-contact couplings would result in more sharply peaked emission 
spectra, due to weaker hybridization and broadening of molecular orbitals.) Since photons 
cannot be emitted at energies higher than eVu as the spectra cut off at this energy. This 
is one of the factors that contributes to the noticeable shift of the emission peak upwards 
in energy with increasing Vu as . Another factor contributing to this blue shift is the cubic 
dependence of the emission rate on the photon energy %u in equation (9). This explains 
the significant photon emission seen in Fig. 2f for Vuas = 6 V at energies well above the 
molecule's nominal HOMO-LUMO gap energy, due to transitions from the higher energy 
tail of O*, and to the lower energy tail of O. 



B. Symmetric molecule-contact couplings; Ep close to HOMO or LUMO. 



In many experimental situations, the HOMO-LUMO gap is not centred at the Fermi 
level of the contacts. Often, either the HOMO or LUMO is close to E F 21 ' 25 ' 28 . In Fig. 3a, O 
has a zero bias energy located just below Ep. In our model, as Vuas increases, O and O* first 
decrease in energy following jiu- Then, at around 3 V, as the energy of O* approaches /is, 
O and O* begin to rise in energy with /is- The opposite situation occurs when O* has a zero 
bias value just above Ep (Fig. 3b). In this case, O and O* first rise in energy with then 
fall with fip, after Vu as = 3 V. This behaviour is nearly identical to recent results obtained 
through self-consistent density-functional calculations for the HOMO and LUMO energies 
of a gold nanowire attached to two gold contacts 29 , suggesting that our model's approach 
to charging can be quite a good approximation for some molecular wires. Fig. 3c shows the 
orbital occupations for the case of O just below Ep in energy at zero bias (same situation as 
Fig. 3a). The total occupation is less than 2 electrons (1.6 electrons) due to the fact that, at 
zero bias, O is very close to Ep, so that its high energy tail is above Ep and is unoccupied. 
As hd decreases, the molecular orbitals drop in energy so that the total molecular charge is 
maintained at its zero bias value. Eventually, as us increases, it approaches the antibonding 
energy, causing O* to start to fill. The orbitals stop dropping in energy, so that any filling 
of O* is accompanied by an emptying of O. Once O is significantly above fip>, its occupation 
becomes constant with further changes in bias, at about 1 electron, as only the source contact 
contributes electrons. Both orbitals now move upwards in energy, causing O* to stop filling 
so that the total charge is still maintained at its zero bias value. The high energy tail of 
O* stays above /is, so that the orbital only fills to about 0.6 electrons. Fig. 3d similarly 
shows the orbital occupations for the case of O* just above Ep at zero bias (the situation in 
Fig. 3b). In this case, the total occupation is more than 2 electrons (2.4 electrons) due to 
the fact that O* is partially occupied at zero bias, because its low energy tail is below Ep. 
O partially empties and O* partially fills as the orbitals enter the electrochemical potential 
window at a bias voltage near 3 V. The occupation of O* becomes about 1 electron, coming 
from the source contact. The occupation of O is greater (1.4 electrons) because the low 
energy tail of O is still below fip,. Fig. 3e shows the total emission rate predicted by our 
model for the case where O is just below Ep at zero bias. A nearly identical emission rate 
is predicted for the case of O* just above Ep. In both cases, emission becomes strong above 
Vbias = 3 V. This is the point at which O and O* enter the electrochemical potential window, 
and therefore are both partially occupied. Emission at high bias is less than half as strong 
as in the case of HOMO and LUMO that are symmetric about E F (Fig. 2), due to the fact 
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that either O* is not as highly occupied (Fig. 3c) or O is not as empty (Fig. 3d), as in Fig. 
2c. Fig. 3f shows the emission spectrum for the case of O just below Ep at zero bias, for 
Hias = 6 V. (A very similar spectrum is predicted for the case of O* just above E F .) The 
spectrum is not blue-shifted as strongly as in Fig. 2f, because the higher energy tail of O* 
remains unoccupied (Fig. 3a) (or the lower energy tail of O remains fully occupied (Fig. 
3b)), resulting in fewer high energy transitions. 

C. Asymmetric molecule-contact couplings; HOMO-LUMO gap centred at Ep. 

By examining the effects of varying /3_i )(J and f3\j, we are also able explore the behavior 
of photon emission from molecular wires with asymmetric molecule-contact couplings, such 
as systems where one of the contacts is an STM tip. In this section, we consider the case 
where the HOMO-LUMO gap of the molecule is centred at the zero bias Fermi level of the 
contacts. First, we consider highly asymmetric couplings: = —0.1 eV, f3\^ = —1.0 eV. 

Fig. 4 shows how O and O* change in energy as a function of bias. The coupling of the 
molecule to the drain is much greater than to the source, so the drain has a much greater 
effect on orbital occupations. As Vu as increases from V, the energy levels of O and O* 
drop at the same rate as jip,- In this way the occupations of O (although a part of its high 
energy tail is above fin) and O* (due to the part of its low energy tail that is below hd) 
change little, and the total charge is maintained at its zero bias value. As O* crosses us in 
energy, its occupation increases very slightly. This causes a slight deviation of the orbital 
energies from a straight line, as O empties slightly by moving closer in energy to fin- The 
occupation of O remains close to 2 electrons throughout, and O* remains almost completely 
unoccupied. In agreement with this, our model predicts that photon emission is extremely 
weak in this situation. If the values of /3-i, a and f3\^ are switched, O and O* rise in energy 
with us instead. The orbital occupations are similar to those in the opposite situation, 
since this time the source has a much greater effect on occupations. So, photon emission 
in this reverse situation is also negligible. This lack of significant photon emission from 
the molecule, caused by the asymmetry of the contact couplings, is a possible explanation 
for the lack of molecular-based photon emission observed in recent STM-HBDC monolayer 
experiments by Hoffmann et al. 37 . 

To understand better the dependence of the photon emission on the asymmetry of the 
couplings, it is useful to consider an intermediate case of asymmetric couplings shown in 
Fig. 5: /3_i )a = —0.6 eV, f3\j, = —1.0 eV. Fig. 5a shows how the model predicts O and 
O* to change in energy as a function of bias. Up to Vbi as = 3 V, O and O* change little 
in energy. Above Vuas = 3 V, both O and O* descend in energy with //p. Fig. 5b shows 
the orbital occupations. For Vu as < 3 V, O* fills significantly as it approaches /is- For this 
reason, at Vuas < 3 V, O and O* do not drop much in energy, as they do in Fig. 4 due to 
the emptying of the high energy tail of O. At Vb ias = 3 V, O* partially fills as it crosses /is, 
but the drain couples more strongly to the molecule. Therefore, O and O* are pulled down 
with //£>, such that O only empties to the extent that O* fills, so that the total charge is 
maintained at its zero bias level. The result is that O empties to about 1.5 electrons, and 
O* fills to 0.5 electrons. Fig. 5c shows photon emission for this case, and other intermediate 
cases. For this case (/3-i, a = —0.6 eV), emission is quite drastically reduced compared to the 
case of symmetric couplings, but is not negligible. For the reverse situation (/5_i j0 = — 1.0 
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eV, Pi t b = —0.6 eV), the reverse arguments apply, and the result is that the emission is 
similar. The degree of contact coupling asymmetry clearly plays a key role in determining 
the strength of photon emission. Our results indicate that in order for an STM experiment 
to detect molecular-based photon emission, the STM tip must not be very weakly coupled 
with the molecule. Since coupling strength varies exponentially with coupling distance, any 
molecular-based photon emission observed in an STM experiment will be very sensitive to 
the tip-sample distance. If the tip is too far from the sample, emission will be negligible. 

In order to facilitate comparison of this prediction with the results of potential experi- 
ments, it is useful to relate contact coupling asymmetry to the amount of current flowing 
through a molecular wire as well, because, unlike coupling strength, current is a directly 
measurable quantity. Fig. 5d shows how the current (obtained from eq. (7)) changes with 
Vbias for various values of /3_i )0 , holding (3\j, fixed at —1.0 eV. As O and O* enter the Fermi 
energy window, these orbitals act as channels for electron transmission, and the current 
increases as a result. Asymmetric contact couplings result in reduced current flow. Notice, 
however, that the current flow is not affected by strongly asymmetric couplings as much as 
is the photon emission rate. Fig. 5e compares how photon emission and current depend on 
contact coupling asymmetry, for Vuas = 6 V. Emission has a power law dependence (roughly 
P-i,a 3 ' 8 ), while the dependence of the current on the asymmetry is clearly much less strong. 
Fig. 5f shows the calculated photon emission vs. the current, for Vu as = 6 V. At weak cur- 
rents, the photon emission rapidly becomes negligible. This is an important consideration 
to keep in mind when attempting to observe molecular photon emission experimentally. 

D. Asymmetric molecule-contact couplings; Ep close to HOMO or LUMO. 

Our results suggest that, for a HOMO-LUMO gap that is symmetric about Ep, asym- 
metric molecule-contact couplings, such as those where one of the contacts is an STM tip, 
result in very low photon emission rates. We consider now strongly asymmetric couplings 
a = -0.1 eV, /3 1>6 = -1.0 eV) for the case of HOMO or LUMO very close to E F at 
zero bias. This coupling configuration is analogous to an STM experiment performed under 
forward bias (the STM tip acts as the (weakly coupled) source, from which electrons flow). 
The opposite configuration (/3_i )0 = —1.0 eV, = —0.1 eV) is analogous to an STM 
experiment performed under reverse bias (the STM tip acts as the (weakly coupled) drain). 
Fig. 6a shows how, for the forward bias situation, O and O* change in energy as a function 
of bias, for the case of O very close to Ep at Vu as = 0. O and O* drop in energy following 
the electrochemical potential fip> of the strongly coupled contact. Even as O* crosses /is, 
the orbital energies continue to decrease almost linearly because, unlike in Fig. 4, O is very 
close to /xd, so that the density of states at //d is high and the weak filling of O* that occurs 
as it crosses /xs can be compensated by a weak emptying of O with only a very small devi- 
ation of O from its linear behavior. Fig. 6b shows the energies of O and O* for the case of 
reversed contact coupling strengths (/3_i j0 = —1.0 eV, f3\^ = —0.1 eV), which is equivalent 
to performing an STM experiment under reverse bias. This time, O and O* rise in energy 
following the electrochemical potential /is- O* never approaches /is, so it remains empty. 
Fig. 6c, d show photon emission for the forward bias situation (Fig. 6a). Biases above 3 V 
result in a small amount of photon emission because, unlike in the case of Fig. 4, O is only 
partially occupied, so some transitions from O* may occur. Emission is still weak (1-2 orders 
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of magnitude less than for symmetric contacts), because 0* becomes only weakly occupied 
when crossing the Fermi level of the weakly coupled contact The emission spectrum, 
shown in Fig. 6d, is peaked at an energy slightly less than 3 eV, because the lower energy tail 
of O is completely filled, and may not receive transitions. In the reversed coupling situation 
(Fig. 6b), photon emission is negligible, since O* never approaches fis and remains empty. 
Now we consider the cases where O* is very close to E F at Vu as = (Fig. 6e,f). Again, the 
energies of O and O* follow the Fermi level of the strongly coupled contact. Therefore, for 
the case of weak source contact coupling (Fig. 6e), there is negligible emission because O 
never crosses \in and remains completely filled. For the case of weak drain contact coupling 
(Fig. 6f), there is (weak) photon emission, because O enters the electrochemical potential 
window and (slightly) empties. (The emission curves for this case are very similar to Fig. 
6c, d.) This property of asymmetrically coupled systems may be used to help distinguish a 
system where the HOMO is close to E F from a system where the LUMO is close to E F , 
with the use of an STM as the weakly coupled contact. Within our model, if the HOMO is 
close to E F , the result is photon emission under high forward bias and negligible emission 
under reverse bias, whereas if the LUMO is close to E F , the result is emission under high 
reverse bias and negligible emission under forward bias. The implications of this may be 
important for molecular electronics research. For example, there is currently disagreement 
over whether E F is located close to the HOMO or the LUMO in the system of BDT attached 
to gold contacts 14 ' 19 ' 21,25 ' 28,30 . The photon emission characteristics of such a system could 
help determine the locations of the HOMO and LUMO, and therefore which type of model 
is more appropriate for such a system. 

IV. CONCLUSIONS 

We have examined theoretically the possibility of photon emission occuring due to molec- 
ular transitions in current-carrying molecular wires. Our results suggest that significant 
photon emission may occur for bias voltages near or above the HOMO-LUMO energy gap 
of the molecule, depending on the specifics of the electronic structure of the molecule and 
contacts and of the couplings between the molecule and the source and drain. The predicted 
photon spectra are peaked near the gap energy 40 . If two molecular orbitals are located in 
the energy window between the electrochemical potentials of the contacts, they will both 
be partially occupied and, as long as transitions between the orbitals are not forbidden, 
transitions from the higher energy orbital to the lower energy orbital should occur, resulting 
in photon emission. We find that varying the strength of the molecule-contact couplings 
changes the emission strength and spectrum. Weak, symmetric couplings result in a sharply 
peaked emission spectrum. Asymmetric couplings result in reduced photon emission. Emis- 
sion from strongly asymmetrically coupled systems (such as some STM-molecular monolayer 
systems) is predicted to be extremely weak if the zero bias value Fermi level E F of the con- 
tacts is not close to the HOMO or LUMO energy of the molecule. If E F is close to the 
HOMO or LUMO energy, the amount of photon emission observed under forward bias and 
under reverse bias, in an STM experiment, may be compared in order to determine which 
of the orbitals is close in energy to E F . 

Whether E F is close to the HOMO or LUMO energy is a difficult unresolved issue in 
current research on molecular wires such as BDT attached to gold contacts; ab initio cal- 
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culations based on different models and approximations 21,25,28 have yielded opposite results 
in this regard. It is believed that the HOMO and LUMO of BDT attached to gold contacts 
are composed of bonding and antibonding orbitals (tt and ir*). The present work suggests 
that, if the contact couplings are symmetric, photon emission will not be drastically differ- 
ent whether Ep is close to the HOMO or to the LUMO: Each orbital will become partially 
occupied as the bias voltage becomes greater than the HOMO-LUMO gap, and, so long as 
7i*— > 7i transitions are allowed, this will result in photon emission. Optical transitions are 
forbidden if (ip f\x.\ipi) = 0. We have used extended Hiickel theory to examine the HOMO 
and LUMO of BDT 42 . The contacts were not included in this calculation. By symmetry 
analysis of the HOMO and LUMO, we found that (HOMO\y\LUMO) ^ 0, where y is in 
the direction of the plane of the molecule, perpendicular to the contacts. Therefore, optical 
transitions should be allowed, and would result in the emission of y-polarized photons. Ex- 
perimental evidence of BDT optical transitions is available in the form of absorption spectra 
for BDT in hexane 43 . Absorption occurs for energies higher than about 4 eV, indicating a 
possible experimental value for BDT's HOMO-LUMO energy gap. Transport measurements 
on Au/BDT/Au molecular wires have been reported at bias voltages as high as approxi- 
mately 5V 1 . Thus it appears to be quite reasonable to undertake an experimental search for 
electroluminescence due to molecular transitions in Au/BDT/Au systems. An experimental 
determination of the HOMO-LUMO gap of BDT in the presence of the Au contacts by this 
means would be an important test of the validity of various theories of molecular wires. 
If one of the Au contacts is an STM tip, electroluminescence measurements may also be 
expected to be helpful in determining the location of the Fermi level of Au relative to the 
HOMO and LUMO states of the molecular wire. 

The study of electroluminescence from molecular wires has only just begun. We hope 
that the theoretical ideas presented here will stimulate further experimental and theoretical 
work on this interesting and potentially important topic. 

This work was supported by NSERC and by the Canadian Institute for Advanced Re- 
search. 
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FIGURES 



FIG. 1. A schematic diagram of the model molecular wire. The source and drain contacts are 
semi-infinite. The two atoms in the center represent the molecule. 

FIG. 2. Symmetrically coupled molecular wires for which at zero bias the HOMO-LUMO gap is 
centred at the Fermi level of the contacts and e C ontads = £a = Q>- a) Source and drain electrochemical 
potentials fi$ and \ib and energies of the bonding (O) and antibonding (O*) orbitals as a function 
of bias voltage, for /3_i i(J = = —1.0 eV. b)Probability for the transmission of an electron 
through the molecular wire as a function of electron energy, for Vu as = and /?_i )0 = = — 1.0 
eV. c)Occupations of the molecular bonding and antibonding orbitals, for /3_i >a = [3\ t b = —1.0 eV, 

d) for /3-i 5 a = Pi,b = —0.2 eV. e)Total integrated photon emission, as a function of bias voltage, for 
(i)/3_i )Cl = 01^ = —0.2 eV and (ii)/3_i i(l = (5\^ = —1.0 eV. f)Emission spectra for various different 
bias voltages (/3-i, a = Pi,b = — 1-0 eV). 

FIG. 3. The case of Ep close to either the HOMO or LUMO at zero bias. /3_i )0 = fi ljb = -1.0 
eV for the entire figure. &)Ep is close to HOMO at zero bias. Energies of the bonding (O) and 
antibonding (O*) orbitals as a function of bias voltage. b)Energies of O and O*, for the case of Ep 
close to LUMO at zero bias, c) Occupations of O and O* , for the situation in (a), d) Occupations of 
O and O*, for the situation in (b). e)Total photon emission rate for the situation in (a) (emission 
rate for situation in (b) is very similar), f) Emission spectrum for the situation in (a) for Vbi as = 6 
V (Similar spectrum for situation in (b)). 

FIG. 4. The energies of O and O* as a function of bias voltage, for the case of a HOMO LUMO 
gap that is centred at the zero bias Fermi level of the contacts, with highly asymmetric contact 
couplings (P- 1>a = -0.1 eV, (3 hb = -1.0 eV). 

FIG. 5. The case of a HOMO-LUMO gap that is centred at the Fermi level of the contacts zero 
bias, for various asymmetric contact couplings. a)/3_i ja = —0.6 eV, (3i jb = —1.0 eV. Energies of 
the bonding (O) and antibonding (O*) orbitals as a function of bias voltage, b) Occupations of O 
and O*, corresponding to the situation in (a). c)Total photon emission rates, for various values of 
P-\,a- fli,b = —1-0 eV in all cases. d)Current for various values of /3-i, a - Pi,b = —1-0 eV in all cases. 

e) Photon emission (solid line) and current (dotted line) as a function of source contact coupling 
/3-i, a , for = —1.0 V, Vbias = 6 V. Vertical scales are arbitrary. f)Total photon emission vs. 
current, for /3 ljb = -1.0 V, Vbias = 6 V. 

FIG. 6. a)The energies of O and O* as a function of bias voltage, for the case of O located 
at Ep at zero bias, with highly asymmetric contact couplings (/3-i, a = —0.1 eV, (3i jb = —1-0 eV.) 
b)Energies of O and O* for reversed couplings (/?-i, a = —1.0 eV, (3i jb = —0.1 eV). c)Total photon 
emission, corresponding to the situation in (a). d)Emission spectrum, corresponding to situation 
in (a), for Vbi as = 6 V. e)Energies of O and O*, for the case of O* located at Ep at zero bias, with 
(3-i : a = —0.1 eV, (3\ t b = —1.0 eV. f)Energies of O and O*, for reversed couplings. 
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